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Page 6, equation (lO) : 
f 

instead of „ 
ft r h 


The second term on the right side should be 


Page 6, equation (ll) : The equation should read as follows: 



Page 7: The second equation following equation (15) should be 

t = -1 instead of r = 0. 

Page 8, equation (18): A bracket should be inserted between J 0 (ctr) 
and C-j-. The bracket should be closed at the end of the equation. 

Page 8, last line: The sentence should read, "A sketch of the stream- 

lines in a radial section is given in figure 5.” 
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TECHNICAL NOTE 4323 


NATURAL CONVECTION INSIDE A FLAT ROTATING CONTAINER 
By Simon Ostrach and Willis H. Braun 


SUMMARY 

Froude number Is found to Be the parameter that determines the ef- 
fect of hody rotations on internal natural convection flows. Explicit 
forms of Froude number are determined for the cases of predominant rota- 
tion and predominant axial force. In a heated cylindrical container no 
significant convective effects are generated by rotations alone. More- 
over, rotation tends to inhibit the flow and heat transfer generated by 
an axial body force in such a configuration. Insertion of radial vanes 
in the container eliminates the detrimental effects of rotation and im- 
proves the heat transfer. 


INTRODUCTION 

A novel application of natural convection to yield a lightweight, 
efficient, and simple heat sink is described in reference 1. The scheme 
described therein, which permits turbine blades to operate successfully 
in a high-temperature gas. Is that a liquid is circulated through cylin- 
drical cavities in idle blades by natural convection generated by centrif- 
ugal force. One might then wonder whether the same Idea might also be 
employed in another difficult problem, namely, that of cooling a high- 
speed vehicle entering the earth's atmosphere. In such a case, two body 
forces are available, the force of retardation by the atmosphere and the 
centrifugal force generated by the rotation of the vehicle about its 
axis. 


Despite the practical demonstration of the basic idea in reference 
1, relatively little information on natural convection exists for inter- 
nal flows and for body forces which vary temporally or spatially. The 
problem of natural convection in a reentering vehicle is further compli- 
cated by the simultaneous action of the two kinds of body forces in dif- 
ferent directions. 

This paper, therefore, presents an analysis of a somewhat idealized 
configuration which, nevertheless, contains some of the complications 
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associated with actual problems. In this way, the primary factors in- 
fluencing the phenomenon are determined, and the relative importance of 
some of the interacting effects can be evaluated. This work thus repre- 
sents an initial study of the use of fluids driven by body forces for 
heat sinks in rotating reentry vehicles. 


DEFINITION OF THE PROBLEM 

A fluid contained in a high-speed vehicle which enters the earth's 
atmosphere is subject to an apparent gravitational force. A represent- 
ative deceleration curve for such a vehicle calculated from the equations 
in reference 2 is presented in figure 1. Clearly, the body force is 
strongly time-dependent. If a representative calculation of the heat 
flux to the vehicle made from reference 2 is superposed on figure 1, it 
can be seen that a lag exists between the heating of the body and the 
appearance of a large retardation force. During the initial period of 
heating there is no driving force for the convective cooling process. 

One wonders whether the heat transfer to the internal fluid might be in- 
creased by rotating the vehicle nose about its axis. Centrifugal force 
would then generate the motion as in the turbine-cooling problem. 

Rotation may be of value even after the decelerating force has be- 
gun to act. The decelerating force ants in the axial direction (see fig. 
2) while the internal fluid near the external stagnation point of the 
blunt nose is heated from below; that is, the heating in the stagnation 
region imposes a negative temperature gradient parallel to the retarda- 
tion force. In such a configuration the fluid remains at rest until a 
critical value of the Rayleigh number (PrGr) is attained (see ref. 3). 
Thus, even for appreciable retardation foraes, the Internal fluid Is in- 
effective as a coolant. Therefore, a large body force (such as a cen- 
trifugal force) transverse to the temperature gradient might Increase the 
effectiveness of the heat sink by starting the motion sooner. Also, the 
rotational effects must be evaluated because rotation of the vehicle may 
be used for aerodynamic stability and, therefore, be inherent to the 
system. 

Specific consideration is given in this report to the flow and heat 
transfer of a fluid subject to an axial body force inside a rotating- 
right circular cylinder of small height (see fig. 3). In this way, the 
rotational effects can be evaluated, and internal conditions in the vi- 
cinity of the external stagnation region of blunt bodies are simulated. 

In order to avoid the complications introduced by unsteady effects, we 
shall assume that the retardation force is constant and concentrate on 
the effects of rotation of the vehicle on the internal flow. Therefore, 
two cases will be considered: (l) where rotation predominates, that is, 

where the internal motion is essentially due to rotation of the container; 
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and (2) where the deceleration force predominates, that is, where a flow 
is generated by the axial "body force with heating irrespective of whether 
the container is rotating. 


ANALYSIS 

The governing equations for the coolant are those expressing the 
conservation of mass, momentum, and energy in cylindrical coordinates 
for compressible, viscous, and heat-conducting fluids subject to a body 
force. With fluid properties constant, these equations are, respectiv ely : 

Continuity: 




3z 


CD 


Radial momentum: 




Azimuthal momentum: 
P 




Axial momentum: 
P 

Energy: 


pc. 


where 


3r^ r Sr r^ 3 9 ^ 3z^ 


( 2 ) 


(3) 
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X = 



_U 

r 


+ 


1 5v 

7 ^. + 


Sz 


All symbols are defined in appendix A. 


Predominant Rotation 

The case of predominant rotation simulates conditions encountered 
by an internal fluid, both during the time shortly after reentry when 
the deceleration force is negligible and in the stagnation region before 
a flow has been established by the instability. 

The state equation 


p = p(T,F) 

for a fluid may be approximated near a reference condition (o) by 

p = p 0 [l - 0(T - Tq) + k(p - p 0 )] 

where p is the coefficient of volume expansion and k is the isother- 
mal compressibility. It can be shown from the conservation of mass In 
the closed vessel to be considered here (fig. 3) that. If the heating 
(temperature difference) on the bottom and the cooling on the top are 
Identical functions of the radius, and if there are no azimuthal or axial 
gradients of press lire, the pressure difference (from the reference state) 
vanishes identically. It Is assumed here that the fluid in the container 
is heated and cooled symmetrically as described and that, even though 
same pressure gradients may exist, they have little effect on the density, 
(if a net heating of the contained fluid is desired, a s mal l bulging of 
the container, properly distributed, will compensate the pressure incre- 
ment.) Then, as for unconfined flows, 

p = p 0 [l - p(T - T 0 )] 

Axlsymmetric flow is assumed in the cylinder of figure 3. The char- 
acteristic velocity of the fluid must be proportional to fiR, and It 
should depend on a representative temperature difference between the 
heated and unheated states. Let Tq be the temperature of the unheated 
state and Tj the temperature at the center of the lower (heated) plate. 
By the condition of symmetric heating and cooling, the temperature at the 
center of the upper plate is Tq - (Tj - Tq) . If the temperature and 
density of the fluid are written as 


T = T 0 + 2(T Z - T 0 )t 
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and 

P = P 0 [l “ 2P(T 2 - T 0 )rj 

it is apparent that the perturbation parameter for the system is 
e = 2p(Tj - Tq). ® ie analysis of appendix B, beginning with a general 

power dependence of the characteristic velocity upon e, shows that only 
a linear dependence is consistent with the dyn ami cs of the problem. 

Thus , the representative velocity is efiR . 

It is now convenient to transform to dimensionless 
rotating system as follows : 

r = Rr, 0=9 + fit, z = hz 

U = eftRu, V = fiR(r + ev) , W = eflR(h/R)w 

P = p c /l 2 R 2 (p 0 + ep) 

The basic equations (l) to (5) to zeroth order in e reduce to 

^Pq _ # 5pq _ fh 

'ST ~ Ti "5T ~ ~ (fiR.) 2 


which imply a solid-body rotation of the fluid when there is no heating 
(e = 0) . To first order in e, that is, the perturbations from solid- 
body rotation due to heating, equations (l) to (5) are 



variables in a 
( 6 ) 
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2 

The ratio Re^/Gr^ is actually a Froude number and represents the 

ratio of inertia to buoyancy forces. In this particular case under con- 
sideration, the Froude number is large; in fact, from the definitions of 
the Reynolds and Grashof numbers based on _G, the reciprocal of the Froude 
number is just e. Therefore, the dominant inertia terms are the buoyancy 
and Coriolis terms, that is, the last two terms on the left side, respec- 
tively, of equation (8) and the last one on the left in equation (9). In 
the energy equation the convection and frictional heating terms are neg- 
ligible with respect to the conduction term so that the heat transfer is 
not affected by the motion and is merely due to conduction. All this im- 
plies, of course, that large velocities cannot be obtained by applying a 
temperature gradient transverse to the centrifugal body force in such a 
rotating configuration. To understand this seemingly unusual result, let 
us solve the equations that result by omitting the negligible terms and 
assuming h/R <<1. Also, the parameter is small for practical 

cases . 


Equation (7) is unchanged, and equations (8) to (ll) become, 
respectively. 


i a 2 u _ 

Re^ Sz 2 " 


+ rr 


2v 


i a 2 v 
Re n 5z 2 


2u 



( 12 ) 


(13) 


( 14 ) 
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T 


dz 2 


= 0 


(15) 


The boundary conditions on the upper and lower surfaces (z = +l/2) 

are 


u = v = w = 0 


At the upper surface, also, 


t = 0 

The the rmal boundaiy conditions at the surface will be discussed 
a posteriori. 

By following Davies (ref. 4) the velocity and temperature distribu- 
tions are written as 

u = a JlO 31 *)^*) 
v = | J 1 (ar)V(z) 

where J 0 (ar) and J^Car) are zero- and first-order Bessel functions, 
respectively, and a is the first root of J-^a) = 0. Qualitatively, 
the assumed radial variations in u, v, and w lead to the loop flow 
shown in figure 4. The viscous boundary conditions at the radial bound- 
ary (end) of the cylinder are satisfied by u and v, but for the w 
component the viscous conditions are replaced by a slip condition. This 
is consistent with the omission of terms of order (h/R) 2 in the differ- 
ential equation. 

The system of equations (7), (12), (13) , (14), and (15) can be com- 
bined to yield an ordinary equation for V: 

(Be a )^V v + 4V’ = -2 
with the boundary conditions 

V(l/2) = V(-l/2) = 0, V* (l/2) = V' (-l/2) = 0, V" (l/2) = 

Solution of this boundary value problem then yields 

u = ~ Ji(a-r) sinh-^/Re^' z cos-^Re^ z - C 2 cosh-^Re^. z 

( 16 ) 


V"(-l/2) = 0 



w = J 0 (ctr)w(z) 
t ~ Ji(ar)T(z) 
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= — J^Car) ^ cosh^/Re^ z sirn/Re^ z + C 2 sinh-^/ReQ z cos^/Re^ z - i zj 

(17) 

v - -4= JoCaxloisinh ^*5 sin ^ + cosh oos 

V Re a ' f 


+ cosh — g — cos — £- 


( -i/ Re o V Re o V Re n 

+ Cg l cosh -1-^ cos -~3i _ slnh sin 


t = — J^arH-z) 

nfr* 


- C 1 ^sir±L- v /Re a z sin^/Re^ z + coshn/Re^ z cosV Re ft 

- C 2 ^coeh-y/Re^ z cos^Re^ z - sinh-^/Re^ z sin-/Re fi z^ (18) 


p = 0 


where 



The dimensionless temperatures of the plates T^r, ± -|j are, from equa- 
tion (19), equal to +(l/ar ),!-]_ (or) . Any other temperature distribution 
could be specified at the lover surface by expanding ar-rfr,- -g) in a 
series in 


The velocity distributions (eqs. (16) to (18)) are presented in fig- 
ure 4. A strong Reynolds number effect is evident. For lov Reynolds 
numbers the viscous effects are present over the entire cross section of 
the cylinder, but at the high Reynolds number a boundary layer appears 
in all three velocity component profiles. . Accompanying the boundary layer 
is a dimi nution of the radial and axial components of velocity in the main 
body of the fluid vhile the azimuthal velocity distribution is nearly 
linear, lagging the rotation of the cylinder in the upper half and leading 
in the lover half. A sketch of the streamlines is given in figure 5. 
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Since the equations of motion have been solved, it is possible to 
explain why, physically, the velocities remain small. The driving force 
rx in equation (12) is opposed by the pressure gradient, Coriolis force 
(-2v), and viscous force. The pressure gradient balances that part of 
the driving force which is independent of z. At low Reynolds numbers 
the remainder of the driving force is opposed mainly by the viscous 
force, but at large Reynolds numbers the Coriolis force grows relatively 
larger and inhibits the flow. Consequently, the heat convection is neg- 
ligible both for large and small Reynolds numbers, and the fluid essen- 
tially rotates as a rigid body. 

Figure 6 is a graphical illustration of the Coriolis force in the 
top of the cylinder. The Coriolis force has negative components in both 
the radial and angular directions. It opposes the azimuthal rotation of 
the fluid by the viscous force and the outward radial flow set up by the 
buoyancy force. 

It may be that flows of a cellular type will be obtained in this 
configuration if h/R is not negligibly small and the condition of axial 
symmetry is relaxed. However, it is doubtful that even then the heat 
cvj transfer would be affected by the fluid motion. Thus, rotating a fluid 

g container like that just considered does not significantly affect the 

heat transfer when the deceleration force is inoperative. In fact, the 
rotation may even be detrimental, for Chandrasekhar (ref. 5) has shown 
that in an unstable configuration the fluid motion is delayed by the ac- 
tion of Coriolis forces. 

Since the Coriolis force due to the circumferential velocity compo- 
nent opposes the radial buoyancy force and thus negates convection effects 
and also may impede the unstable motion, it appears that significant con- 
vection could be obtained if the Coriolis forces were made negligible 
compared to the centrifugal force. The insertion of radial vanes appro- 
priately spaced would arrange the circumstances by reducing the azimuthal 
velocity component (and, hence, the corresponding Coriolis _ term) . The 
buoyancy effects due to centrifugal force would, in this way, be unop- 
posed (excepting viscous forces) and could lead to large heat convection 
as in the turbine blades investigated by Schmidt (ref. l) . 

Consider a sector between two vanes in such a divided flat cylinder 
(fig. 3(b)). It is anticipated that, in contrast to the previous con- 
figuration, the inertia forces will be of the same order of magnitude as 
the buoyancy forces. It follows that the representative velocity must 
be -i/eftR rather than eftR. Then, upon the equations of motion (eqs. 

(l) to (5)) the following transformation to a dimensionless, rotating 
system is performed: 


r = Rr 
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6 = (&p)y + at 

~z = hz 
U = -^sORu 
V = ORr + ViimC^p)v 
V = -/eClhw 

The new angul ar variable y takes values from 0 to 1. Again, the 
state variables are given by 

P = p 0 a 2 R 2 (p 0 + ep) 

P = P 0 (l - €t) 

T - T 0 = 2t(T z - T 0 ) 

The transformed equations, with the small compressibility and dissipation 
terms omitted for simplicity, are 

No heating: 

POr = r 

p_ = 0 

■^Oy 

fh 

P ° z ” " (OR) 2 

Continuity: 

(ur) r + v y + (rw) z = 0 

Momentum: 

uUy + ~ Uy + wUj, - G ^) 2 + rT - 
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If the sector Is made slender "by setting 

| « 1, Acp « 1 

the momentum and energy equations reduce to 



In these equations the inertia and heat convection terms are sig- 
nificant in contrast to the previous case of the complete cyl inde r. The 
presence of the nev parameter £ep controls the Coriolis force and makes 
this possible. The magnitude of the Coriolis terms is dete rmin ed by the 



12 


NACA TN 4323 


square root of a Proude number Acp/ ~/e } and the ratio h/ (R A?) controls 
a portion of the viscous force. Making A*p very small reduces the 
Coriolis force hy preventing the acceleration of the fluid to appreciable 
azimuthal velocities relative to the rotating system. At the same time, 
however, the ratio h/ (R Acp) may become large, which indicates the flow 
has become constricted and causes large viscous retardation. Probably 
an optimum geometrical arrangement would be __ 

| < Acp < Ve 

The previous equations resemble those for natural convection over 
a heated vertical flat plate. However, the driving force here 1 b the 
centrifugal force, which can be many times larger than the earth's grav- 
itational force. Hence, the convective heat transfer in such a rotating 
configuration can be much larger than that due to gravitational effects, 
as was demonstrated in the Schmidt turbine (ref. 1). 

It can be seen from the equations that the Coriolis and viscous 
forces in a rotating system will always have their adverse effect. All 
this means, however, is that convective heat transfer produced by cen- 
trifugal forces in a rotating sector can never be as large as that over 
a flat plate in a gravitational field of corresponding magnitude . 


Predominant Deceleration Force 

The effects of rotation on a flow generated by the deceleration 
force must still be evaluated because the reentry vehicle may be rotating 
about its axis for aerodynamic stability reasons. The deceleration force 
will establish flows in a flat cylindrical configuration (fig. 3) after 
the critical Rayleigh number is surpassed. Also, if the retardation force 
is tr an sverse to a temperature gradient, as it would be, for example, in 
a curved container like a spherical shell, it would Immediately generate 
a conventional natural convection flow. Let us, therefore, analyze the 
case where the flow is essentially due to the retardation force. 

For the inertia and buoyancy forces to be of equal order of magni- 
tude, the characteristic velocity must be V efh. Accordingly, make the 
following transformation from the inertial coordinates to dimensionless 
rotating coordinates : 

U = -y/efhu V = -/efhv + DRr 

r = Rr 9 = cp + fit 

P = P 0 (! - 6T) 



P = p 0 fh(po + ep) 


t = (T - T 0 )/2(T Z “ T 0> 
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In general , u, v, and w are functions of the angular velocity as 
well as the coordinates. Note that, if H is zero, the flow is due en- 
tirely to the action of the axial body force and the heating. On the 
other hand, for no heating the flow is again a solid-body rotation. This 
is clear from the equations since, for e = 0, equations (2), (3), and 
(4) become, respectively, 


_ Cor ) 2 
Sr fh T } 




When e ^ 0, and again when small compressibility and dissipation terms 
are omitted, the equations of motion become, to lowest order in e, 


S 

37 


(ru) + 


Sv + s 

"Sep 35 


(rw) 


0 


( 20 ) 


Su , v chi 

u 37 + 73$ 


+ w 



2<nv 



1 R d 2 u /h\ 2 /cl 2 u 1 Su _1_ S 2 u u 2 Sv\ 

-y/oir ^ Sz 2 ^Sr 2 r ^ r 2 Sep 2 r 2 r 2 


( 21 ) 
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In these equations. 


CD = 


OR 


-/Gr = 


The parameter CD is proportional to the Froude number. In the case 
being considered here cd assumes small values by definition because the 
numerator depends on the rotation and the denominator depends on the 
(predominant) retardation force, (in the previous case, the inertia and 
buoyancy effects both depended only on the rotation with the consequence 
that the Froude number was independent of rotation.) Therefore, with 
cd "S, 1, the inertia terms in equations (2l) and (22) having cd as a co- 
efficient are, at most, of the same order as the other like terms. The 
equations look like the conventional boundary- layer equations, and the 
inertia and convection effects are important. The heat transfer will be 
greater than that due to conduction alone. 


An examination of the foregoing equations shows that a steady cellu- 
lar flow, as occurs with heating a horizontal surface from below (ref. 3), 
is affected by a small rotation to a different order of magnitude than an 
axisymmetric flow. Consider first the cellular flow. With fi = 0, the 
three velocity components of each cell are of the same order of magnitude. 
The first momentum equation (eq. (21)) Bhows that the rotation will affect 
this flow to first order in cd. The preceding analysis of the flat cyl- 
inder with rotation only and the analysis of the unstable static config- 
uration (ref. 3) suggest that the effect will be adverse. Thus, it is 
expected that, in cellular flow generated by an axial body force, rota- 
tion will reduce the heat transfer. 


Second, consider a steady axisymmetric flow with zero azimuthal ve- 
locity component. When a small rotation is imposed, the azimuthal veloc- 
ity is properly represented by 


v = (DV^ 

Since the property of axisymmetry is retained (5/dq> = 0) , it is evident 
that the equations are affected (through the Coriolis and centrifugal 
force terms) only to order (S? . Consequently, the rotation will have 
negligible effect upon this type of natural convection flow. 


CONCLUSIONS 

The analysis presented herein indicates that the Froude number is 
the parameter that determines the relative effect -of body rotation on 
internal natural convection flows. In flows with large Froude number 
(e.g., motions induced solely by rotation of a heated cylinder), the 
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action of the Coriolis force prevents appreciable convection of heat. 
Although the Coriolis force inhibits convective heat transfer somewhat 
even for Froude numbers of unit order of magnitude (rotating sector), 
considerable convection is obtained because the driving force can be 
greatly increased. Flows that have been set up by a strong external 
body force have a small Froude number and are not greatly affected by 
rotation of the container, especially if they have axial symmetry. 

A practical conclusion to be drawn from these results concerns the 
use of large centrifugal forces to drive natural convection cooling. 

Such a scheme will be successful only if the Coriolis force is reduced 
by limi ting the azimuthal extent of the cooling passage. This, in turn, 
increases the viscous forces. Nevertheless, such cooling arrangements 
can greatly affect heat- transfer rates, yielding rates many times higher 
than those using the gravitational field, as shown by Schmidt's tests on 
turbine blades. 


Lewis Flight Propulsion Laboratory 

National Advisory C ommi ttee for Aeronautics 
Cleveland, Ohio, June 17, 1958 
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APPENDIX A 
SYMBOLS 

See each section for specific definitions of the nondimens ional 
quantities. ■ ■ 

A acceleration of vehicle 

Cp specific heat at constant pressure 

c v specific . heat at constant volume 

f negative of z- component of "body force per unit mass 

Gr Grashof number, 2(3 - Tg)fh^/v 2 

Gr^ Grashof number due to rotation, (T^ - Tq)/v^ 

H heat transferred per unit area per unit time 

h height 

frictional heating parameter, ^(^®) 2 /Cp 
k thermal conductivity coefficient 

P pressure 

Pr Prandtl number, Cpp/k 

p nondimens ional pressure 

R radius 

Re Reynolds number due to rotation, fth 2 /y 

*« 

r nondimens ional radial coordinate 

r radial coordinate 

T temperature 

t time 

U,V,W velocity components in radial, azimuthal, and axial directions, 
respectively 
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£ 


U ,v,w 


u,v,v 

y 

z 

z 

p 

r 

€ 

e 

n 

V 

p 

T 

4 

<P 

ft 


dimensionless velocity profiles 
nondimens ional velocity components 
dimensionless angular coordinate in sector 
nondimens ional axial coordinate 
axial coor din ate 

fluid volumetric expansion coefficient, p[- 


5 (1 /P) 

3f 


P 


ratio of specific heats 
constant, 2p(l^ - Tq) 

angular coordinate in inertial coor dina te system 
absolute viscosity coefficient 
kinematic viscosity coefficient 
density 

dimensionless temperature difference 
dissipation function 

angular coordinate in rotating coordinate system 
angular velocity 


<d ftR/-/efh 

Subscripts : 

l center of lover surface 

max maximum 


0 


unheated state 
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APPENDIX B 

ORDER OF MAGNITUDE OF VELOCITY COMPONENTS 

Instead of the relations (6), let a more general dependence of the 
velocity components on e he employed: 

U = flRe m u V = OR(r + e n v) W = QRe^v ' 

r = Rr z = hz P = p^2^R 2 (pq + e^p) (Bl) 

p = p Q (l -ex) t = (T - T 0 )/2(T l - Tq) m,n,l > 0 ! 

Equations (Bl) are appropriate for a flow which increases with external 
heating and reduces to solid-body rotation when there is no heating. 
Substituting into the governing equations (l) to (5) , letting h/R -* 0, 
and neglecting compressibility and dissipation terms yield 

2 

e^Cuu,. + wc^) - e 2n + erx - e n 2v = -e l p r + e^e" 1 ^^ 
e m+n (uv r + wv z + uv/r) + e m 2u = e^Re^v^ 

P z = 0 

In the first of equations (B2) , no term can be of greater order of mag- 
nitude than the driving force err. Hence, m = n = l =1. Then equations 
(Bl) become identical to (6) and equations (B2) become 

eCuUp + wu z ) + rx - 2v = -p r + \ 

e(uv r + wv z + uv/r) + 2u = Re" 1 v zz > (B3) 

P z = 0 / 

If any value other than 1 is chosen for the indices m and n, the 
Coriolis term in either one or the other of equations (B3) will dominate, 
and a null solution will result. The energy equation becomes 

(e/r)(uT r + wx z ) = (PrRe fi )" 1 x zz 



(B4) 
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From the definitions 


fih 


cflV 


Re n - v 3 G:t h ” v 2 


it follows that 


e = 


Sr n 

Re n 


This accounts for the form of equations (8) to (ll) . 
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Figure 1. - Heat flux and deceleration as function 
of altitude for a representative reentry-type 
vehicle . 
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Figure 2. - Cooling shield in a decelerating nose cone. 
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(b) Azimuthal component. 

Figure 4. - Continued. Velocity profiles in rotating cylinder. 
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(c) Axial component. 

Figure 4. - Concluded. Velocity profiles in rotating cylinder. 







